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С использованием неэмпирического квантово-химического метода О/Т (теория функционала плотности) изучены основные 
энергетические эффекты при образовании комплексов органических веществ с белковыми молекулами, которые определяют 
термодинамику этого процесса. Комплексообразование рассматривается как замещение воды в сольватированных молекулах 
белка и органического соединения. Показана определяющая роль вкладов сильных водородных связей и гидрофобного эффек- 
та в прочности образованных комплексов. 


Оценка фармакологической и биологической 
активности веществ имеет огромное практическое 
значение. Целью нашей работы было количествен- 
ное описание энергетических эффектов образова- 
ния комплексов лекарственных веществ с молекула- 
ми белка в условиях взаимодействия в водной среде. 

Термодинамическое описание позволяет рассма- 
тривать образование таких комплексов как динамиче- 
ское равновесие различных энергетических состоя- 
ний. Расчет энергий этих состояний возможен с при- 
менением современных методов квантовой химии. 

Материалы и методы 

Неэмпирические квантово-химические расче- 
ты ОРТ (теория функционала плотности) выпол- 


нялись с помощью стандартного пакета программ 
САМЕ88 методом ВЗРУР/6-3110* с обменным 
функционалом Беке ВЗ и корреляционным функ- 
ционалом Ли, Янга и Пара (БУР) [1]. Геометрии 
рассчитанных молекул были полностью оптимизи- 
рованы, отсутствие мнимых частот колебаний под- 
тверждало их стационарный характер. 

Результаты и обсуждение 

Константу прочности комплекса лекарственно- 
го вещества с молекулами белка в водной среде 
можно оценить по термодинамическому уравне- 
нию: 

-А О -А Н Д5 
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Энтропийная составляющая имеет величину 
Т- ДХ<4 кДж/моль [2] . 

Процесс образования комплекса мы рассматри- 
ваем как динамическое равновесие двух энергети- 
ческих состояний в условиях организма (іп ѵіѵо) 
при 310 К (37 °С). Исходное состояние - гидрати- 
рованные молекулы лекарственного вещества (Л В) 
и белка. Второе состояние - образовавшийся ком- 
плекс молекулы ЛВ с белком и взаимодействие вы- 
свободившихся молекул воды между собой. 

Энтальпийную составляющую ур. (1) можно 
рассчитать по формуле: 

&Н = {&Н ЛВ БС + а Н в в ) — (д Н лв в + а Н БС В ), (2) 

где: А Н ЛВ _ БС - энергия взаимодействия активной 
группы лекарственного вещества с группами моле- 
кулы белка; А Н в _ в - энергия взаимодействия моле- 
кул воды между собой; А Н лв _ в - энергия гидратации 
активной группы лекарственного вещества; А Н БС ^ В - 
энергия гидратации активной группы белка. 

Межмолекулярное взаимодействие описывает- 
ся многими составляющими. Слабые Ван-дер-Ва- 
альсовы взаимодействия имеют энергию менее 
0,4 кДж/моль, а с учетом того, что в исходном состоя- 
нии молекулы сольватированы, т.е. уже взаимодей- 
ствуют с другими молекулами, разница этих энергий 
имеет на порядок меньше значения и близка к нулю. 

Особое место в межмолекулярных взаимодей- 
ствиях занимает водородная связь, которая может 
иметь энергию от 1 до 200 кДж/моль, сравнимую с 
энергиями ковалентных химических связей. Раз- 
личают слабую водородную связь с энергией до 
20 кДж/моль, среднюю - 20. ..60 кДж/моль, обра- 
зующуюся при взаимодействии нейтральных моле- 
кул, и сильную водородную связь от 60 до 
200 кДж/моль, которая наблюдается при взаимо- 
действии молекул с заряженными ионами [3]. 

Был проведен расчет 40 различных органиче- 
ских и неорганических молекул, включая 20 ами- 
нокислот, встречающихся в организме. Также нами 
был проведен расчет димеров и тримеров аминоки- 
слот, который показал, что при полимеризации эф- 
фективные заряды на активных центрах аминоки- 
слотных остатков практически не изменяются при 
наращивании цепи. Это позволяет моделировать 
взаимодействие с белком, заменяя его взаимодей- 
ствием с аминокислотными остатками. Для всех 
рассмотренных молекул была рассчитана энергия 
взаимодействия сильно заряженных центров с мо- 
лекулами воды. 

Затем были рассчитаны энергии замещения для 
процесса комплексообразования по ур. (2). В ходе 
этих процессов число водородных связей не изме- 
нялось. 

Энергетические эффекты межмолекулярного 
взаимодействия в корне изменяются при учете то- 
го, что взаимодействие происходит в водной среде. 
Это иллюстрируют следующие примеры: 


Таблица 1. Энергетический эффект взаимодействия карбо- 
нильной группы с гидроксильной группой с обра- 
зованием водородной связи в водной среде 


Взаимодействие 

Энергия, кДж/моль 

СНз-СООН ... НОСИ - СНз 

-59,85 

Н 2 0 ... Н 2 0 

-43,47 

СНз-СООН ... Н 2 0 

-63,26 

сн 3 -сн 2 он ... Н 2 0 

-43,27 

Энергетический эффект замещения воды 
СНзСООН ... Н 2 0+СН 3 СН 2 0Н ... Н 2 0= 
=СН 3 СООН ... Н0СН 2 СН 3 +Н 2 0 ... Н 2 0 

+3,21 


В данном случае происходит взаимодействие 
карбонильной группы с гидроксильной, если взаи- 
модействие в вакууме, то энергетический эффект 
составил бы -59,85 кДж/моль, в водной же среде 
эти группы уже связаны с молекулами воды с энер- 
гией -63,26 и -43,27 кДж/моль. В результате энер- 
гетический выигрыш имеет величину +3,21 кДж, 
т.е. с учетом энтропийной составляющей, близок к 
нулю. Т.е. введение слабых групп в принципе не 
способно повысить аффинитет Л В. 


Таблица 2. Энергетические эффекты взаимодействия ме- 
тильных радикалов между собой и с молекулами 
воды 


Взаимодействие 

Энергия, кДж/моль 

СН 4 ... СН 4 

=-0,1 

Н 2 0 ... Н 2 0 

-43,27 

СН 4 ... Н 2 0 

=—0,1 


Второй пример иллюстрирует гидрофобный 
эффект. В данном случае выигрыш обусловлен вы- 
сокой энергией взаимодействия молекул между со- 
бой, что и делает этот процесс замещения энерге- 
тически выгодным. Гидрофобный эффект пропор- 
ционален площади гидрофобного сегмента и со- 
ставляет величину -10 кДж/моль-нм 2 , для металь- 
ной группы около -2 кДж/моль, бензольное коль- 
цо дает эффект -11 кДж/моль. В вакууме такое 
взаимодействие имело бы энергию порядка 
-0,1 кДж/моль [4]. 


Таблица 3. Энергетический эффект взаимодействия карбок- 
сильной группы с карбоксильным ионом с обра- 
зованием водородной связи в водной среде 


Взаимодействие 

Энергия, кДж/моль 

СНз-СООН ... ООС- СНз 

-139,24 

Н 2 0 ... Н 2 0 

-43,27 

СНз-СООН ... Н 2 0 

-63,26 

сн 3 -соо~ ... Н 2 0 

-84,37 

Энергетический эффект замещения воды 
СНзСООН ... Н 2 0+СН 3 С00~ ... Н 2 0= 
=СН 3 СООН ... ~00ССНз+Н 2 0 ... Н 2 0 

-34,88 


Третий пример иллюстрирует энергетический 
выигрыш при образовании сильных водородных 
связей. Хотя эти группы имеют большую энергию 
взаимодействия с молекулами воды, тем не менее, 
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энергетический выигрыш процесса замещения со- 
ставил -34,88 кДж/моль. 

Процесс образования прочного комплекса 
можно рассматривать не только как замещение мо- 
лекулами ЛВ воды на аминокислотах, но и как за- 
мещение аминокислотными остатками воды в 
комплексе ЛВ-вода. Чтобы активный центр моле- 
кулы ЛВ заместил молекулу воды, связанную с мо- 
лекулой белка, необходимо, чтобы он был спосо- 
бен образовывать более сильную водородную 
связь, чем молекула воды. И наоборот, чтобы ак- 
тивный центр аминокислотного остатка заместил 
молекулу воды, связанную с веществом, он должен 
образовывать более сильную водородную связь. 
Необходимо, чтобы группы в молекулах лекар- 
ственного соединения и аминокислотах были од- 
новременно сильнее воды, тогда такое замещение 
может усиливать активность препаратов. Однако 
не во всех системах, соответствующих этому пра- 
вилу, наблюдается энергетический выигрыш, то 
есть это необходимое но не обязательное условие. 
Все остальные активные центры как молекул дей- 
ствующего вещества так и молекул белка, образую- 
щие менее сильные водородные связи, в принципе 
не способны дать энергетического выигрыша при 
образовании комплекса в водной среде, то есть уве- 
личить его прочность. 

Этот факт позволяет исключить множество ва- 
риантов взаимодействия, не дающих увеличения 
прочности комплекса в водной среде. В молекуле 
ЛВ следует рассматривать только несколько цен- 
тров, которые дают энергетический выигрыш. Сре- 
ди 20 аминокислот, встречающихся в организме, 
способны образовывать более сильные, чем моле- 
кулы воды, водородные связи только аргинин 
(Ащ), лизин (Бу$), аспарагиновая (Акр) и глутами- 
новая (Сіи) кислоты. При наличии сопряженных и 
ароматических систем в исследуемых молекулах 
водородная связь может иметь аномально высокие 
значения. Среди аминокислот, содержащих сопря- 
женные системы, энергетический выигрыш в вод- 
ной среде могут иметь тирозин (Туг), триптофан 
(Тгр) и заряженная форма гистидина (Шк). 

Был произведен квантово-химический расчет 
модельных систем. В табл. 4 приведены величины 
изменения внутренней энергии в процессах ком- 
плексообразования между органическими молеку- 
лами в водной среде. 

Следует отметить, что энергетический вы- 
игрыш наблюдается только при взаимодействии 
органических анионов с нейтральными молекула- 
ми. Особо следует отметить образование прочных 
водородных связей в ароматических и сопряжен- 
ных системах. Катионные формы молекул хотя и 
образуют прочные водородные связи, однако энер- 
гетических выигрышей не наблюдается, так как ис- 
ходное состояние при их взаимодействии с водой 
велико по энергии. В системах, содержащих только 
нейтральные молекулы, энергетический выигрыш 


незначительный, либо отсутствует вовсе. Все эти 
факты позволяют говорить о комплиментарности 
при взаимодействии органических молекул с бел- 
ком. Только определенные сочетания активных 
групп органической молекулы при их взаимодей- 
ствии с активными группами аминокислотных ос- 
татков способны обеспечить комплексообразова- 
ние. Только некоторые группы, входящие в состав 
молекул действующего вещества, способны повы- 
шать энергию процесса комплексообразования, 
т.е. увеличивать константу прочности. 


Таблица 4. Энергетические эффекты замещения в водной 
среде 


Система 

АЕі, 

кДж/моль 

Система 

АЕі , 

кДж/моль 

НСООН+НСОСГ 

-31,19 

Фенол + НСОСГ 

-48,72 

СНзСООН+СНзСОСг 

-34,88 

Фенол + СНзСОСГ 

-47,83 

СН 3 СООН+СН 3 СООН 

34,87 

Фенол + Фуран 

12,61 

С 2 Н 5 ОН+СН 3 СООН 

3,21 

Фенол + С 2 Н 5 МНз + 

33,03 

С 2 Н 50 Н+СНЗСОСГ 

- 0,20 

Пиррол + СНзСОСГ 

-37,76 

СН 30 СНЗ+СНЗС 00 Н 

-0,82 

С 6 Н 5 СООН + СН 3 СОО- 

-40,11 

С 2 Н 5 МН 2 +СНзСОСг 

28,01 

Фуран + НСООН 

19,05 

С 2 Н 5 МН 2 +СН 3 СООН 

13,92 

Пиррол + НСОО~ 

-37,12 

Пиридин+СНзСООН 

-8,32 

Пиррол + Фуран 

8,52 

Анилин+СНзСОСГ 

-31,24 

Фуран + С 2 Н 5 МН 3 + 

22,31 

Анилин+СНзСООН 

15,82 

Индол + СН 3 СОО~ 

-53,45 


Вторым эффективным способом увеличения 
прочности комплекса является наращивание ги- 
дрофобных участков молекул лекарственных ве- 
ществ в соответствии с гидрофобными участками 
молекул белка. 

Эти примеры показывают, что энергетические 
эффекты процесса образования комплекса в корне 
изменяются при учете того, что взаимодействие 
происходит в водной среде. 

Этот факт позволяет более рационально подхо- 
дить к проектированию лекарственных веществ, в 
значительной мере сокращая затраты времени и 
средств для предсказания активности. Вычисляя 
энергии гидрофобного взаимодействия и рассчи- 
тывая энергетические эффекты замещения моле- 
кул воды с образованием сильных водородных свя- 
зей, мы способны рассчитать максимально воз- 
можный аффинитет данного лекарственного веще- 
ства к идеально подходящему белку, что в свою 
очередь позволяет определить максимально воз- 
можную константу прочности комплекса лекар- 
ственного вещества с белком. Величина макси- 
мального аффинитета ЛВ к белковому субстрату 
служит количественной характеристикой, связан- 
ной с его фармакологической активностью, она 
позволяет оценить необходимые количества ве- 
ществ, чтобы наблюдать фармакологический эф- 
фект, а также токсические дозы. 

Например, произведем оценку максимально 
возможного аффинитета для молекулы фуросеми- 
да, рисунок. 
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Рисунок. Молекула фуросемида 


В этой молекуле имеется несколько групп, спо- 
собных образовать сильные водородные связи с 
белковыми мишенями: сульфаниламидная, кар- 
боксильная и аминогруппа. Однако, необходимо 
учесть, что аминогруппа образует внутримолеку- 
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